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ABSTRACT 

We use the C24668, Fe4383, H7A and HSa spectral absorption line indices, together 
with U- and V-band photometry of 101 galaxies in the Coma cluster, to investigate 
how mean age and metal abundance correlate with galaxy luminosity. In particular, we 
use the line index measurements to address the origin of the colour-magnitude relation 
(CMR). We find that the CMR in Coma is driven primarily by a luminosity-metallicity 
correlation. We additionally show evidence for a relation between age and luminos- 
ity, in the direction pre dicted by the semi-analytic hierarchical clustering models of 
Kauffmann & Chariot ( 1998| ), but this is only present in the C24668 index models, 
and could be an effect of the lack of non solar abundance ratios in the Worthey models 
used. 

By comparing deviations from the CMR, with deviations in absorption index from 
analogous 'index-magnitude' relations, we find that colour deviations bluewards of the 
mean relation are strongly correlated with the hydrogen Balmer line series absorption. 
We show that the properties of these blue galaxies are consistent with the presence of 
a young stellar population in the galaxies, rather than with a reduced metallicity. 

Key words: galaxies: abundances - galaxies:clusters:individual:Coma - galax- 
ies: formation - galaxies: elliptical and lenticular - galaxies :starburst 



1 INTRODUCTION 

It has been long understood that the colours of early-type 
galaxies are governed primarily by the effects of age and 
metallicity, which, when increased, cause the spectral en- 
ergy distr i bution (SEP) to become redder (|Rcnzim fc Bnz- 
zoni 198q; |Buzzoni et al. 1992|; Buzzoni et al. 199S ; Worthey 



1994| ; pharlot fc Silk 1994^ In his '3/2' law (Alog(t) 
|A[Fe/H]), Worthey (1994) demonstrated how age and 



metallicity have a degenerate effect on galaxy colours. The 
colour-magnitude relation (CMR), in which the colours of 
early-type galaxies become progressively redder with in- 
creasing luminosity, and hence i ncreasing mass of the galaxy 
(Visvanathan & Sandage 1977) is seen in the cores of rich 



in which more massive galaxies can retain their supernova 
ejecta for longer than can smaller galaxies, thus being able 
to process a larger fraction of their gas before it is expelled 
from the galaxy. Given the degeneracy between the metallic- 
ity and age of a stellar population in its broad band colours, 
it is also possible that the CMR is an age driven sequence, 
with the smaller galaxies being bluer due to a younger mean 
age of their stellar populations. An age dependant CMR 
however ne ither preserves its slope, nor its magnitud e range 
with time (Kodama 1997; Kodama & Arimoto 1997). Stud- 



clustersr in grotTE ggd ovt j 11 ^J*is tn be present in field el- 



ies of the CMR in high redshift clusters find a ridge-line slope 
comparable to that of the local clusters; furthermore, there 
is no sign of a change in th e range of magni t ude over which 
the C MR may be traced ( Ellis et al. 1997 ; Stanford et al 



lipticals-( Larson et al. 190C ). TiadiliuiiaUy, the si out, 1 seen m 



t he CMR has been attribut ed to a mass-metallicity sequence 
(Dressier 1984; Vader 1986), with the massive galaxies being 
more metal rich, and thus redder, than the less massive ones. 
This tendency can naturally be explained by a supernova - 
driven wind model (Larson 1971; Arimoto & Yoshii 1987), 



199S ) . This makes an age dependant CMR in clusters highly 
unlikely. Given this metallicity driven interpretation of the 
CMR, its low levels of scatter in cluster cores implies that 
the gala xies are made up f rom uniformly old stellar popu- 

Even small 



lations ( Bower et al. 1992] ; 



Bower et al. 



1998] ). 



Email: alc@star.sr.bham.ac.uk 



variations in the ages of the galaxies would lead to unaccept- 
able levels of scatter in young stellar populations, whereas 
old stellar populations have a much smaller age dependency 
in their colours. 
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On closer inspection however, this picture runs into 
some problems. At high redshifts, the fr action of blue galax- 
ies i n many of the clusters in creases (Butcher fc Oemler 



1976 



I utcher fc Oemler 1984). Couch & Sharpies (1987 



spectroscopically investigated the galaxies in some of these 
'Butcher-Oemler' clusters. They split their sample into blue 
and red galaxies depending on whether they lie on the clus- 
ter CMR or not, and associate the increased fraction of blue 
galaxies with either ongoing or recent star formation. They 
also found that 11 out of their 73 red galaxies showed en- 
hanced Balmer absorption lines, indicative of recent bursts 
of star formation, although it is possible that these red H<5 
strong galaxies can be for med by the truncation of st ar for- 
mation in spiral galaxies ( Barbaro fc Poggianti 1997 ). 



Spectroscopic studi es of e arly-type galaxies in local rich 
cluste rs (Caldwell et al. (1993), hereafter C93; Caldwell et al. 
(1996)) have shown that these low z clusters also contain a 
population of galaxies with abnormal spectra. They exhibit 
enhanced Balmer absorption lines, indicative of recent star 
formation, but with too weak an [O II] line to be classified 
simply as spiral galaxies. They note the similarity of some 
of these spe ctra w ith the 'red' H5-strong galaxies of Couch 
& Sharpies ( 1987 ), and find that of the galaxies associated 
with the dynamically separate group of galax ies centered 



on NGC4839 in the SW of the Coma clu ster (Baier 1984 



Escalera et al. 1992 ; Colless fc Dunn 1996 ), about 1/3 have 
abnormal spectra. 

These spectroscopic studies of local and distant clusters 
have shown there to be a population of galaxies which, al- 
though they appear photometrically old, have been forming 
significant quantities of stars in their recent past. 

In an effort to disentangle age a nd metallicity effects 
on stellar populations, Worthey (1994) (also Worthey & Ot- 
taviani, 1997) developed a series of stellar population spec- 
tral synthesis models, to show the dependence on age and 
metallicity of spectral features such as the Balmer lines and 
various metal lines. In a study of a magnitude limited sam- 
ple (Mb < —17) of ear ly-typ e galaxies in the Fornax cluster, 
Kuntschner & Davies (1995 ) used these models to show that, 
as predicted by the high-z studies, the elliptical galaxies are 
indeed uniformly old, and span a range in metallicities. They 
also showed that the majority of the (luminosity weighted) 
young galaxies in the cluster were low luminosity lenticular 
systems. However, earlie r studies of m ainly field early-type 
galaxies ( ponzalez 1993 ; Trager 1997) show just the oppo- 
site, with these galaxies seeming to have a uniform metal- 
licity and a range of ages, yet still forming a well-defined 
CMR. This points to the possibility of an environmental 
dependence of the CMR, or maybe two completely differ- 
ent mechanisms operating in cluster and field environments, 
conspiring to produce similar CMRs. 

In this paper we will look at spectral line indices and de- 
viations fr om th e CMR in the Coma cluster, using both the 
Worthey ( |l994| ) and Worthey & Ottaviani ( |l997p spectro- 
photometric evolution models. We use the high precision 
photo metry of the Coma cluster prese nted in Terlevich et al., 
([I999|) (hereafter T99, also Terlevich (|1998;)) and the spectra 



of Coma cluster galaxies from C93. The selection of galaxies 
and the spectroscopic and photometric measurements are 
outlined in §2. Our principle goal is to investigate whether 
the CMR in Coma is principally driven by age or metallicity 
and to determine the cause of deviations from the CMR in 



individual galaxies. Our analysis is presented in §3, followed 
by a summary of our conclusions in §4. 



2 THE DATA 

The overall aim of this paper requires us to combine accu- 
rate photometry with precise line-index measurements. The 
aperture photometry of T99 covers an area of the Coma 
cluster that is well matched to the spectroscopic survey of 
C93. As a result, we are able to compare the spectroscopic 
and broad-band colours of 101 galaxies drawn from within 
50' of the cluster center^]. Below we describe the selection 
of the spectroscopic catalogue, the spectroscopic determina- 
tion of line strength indices, and the matching photometric 
measurements. 



2.1 Sample selection 

C93 selected their sample from the exte nsive galaxy cata- 
logue of Godwin, Metcalfe & Peach ( |l983| ) (hereafter CMP). 
The parent catalogue is drawn from a 2.6 deg 2 field, centered 
on the Coma cluster, and is considered complete to B = 20, 
corresponding to Mb < —14.20 . C93 used two criteria to se- 
lect early-type galaxies from the CMP sample. Firstly they 
fitted the GMP (B, B — V) |§] colour magnitude relation and 
selected all galaxies whose colours lay within ±0.15 mag 
of the fit. Secondly, they obtained independent morpholo- 
gies using KPNO 4m and Palomar Sky Survey plate s. Us- 
ing these morphologies, and those of Dressier ( |i98oD the y 
removed morphologically late-type systems from the list, 
and added some (~ 25) morphologically early-type galax- 
ies which had failed to be included after the colour cut. 
Finally they rejected all galaxies brighter than B = 14.3, 
as high velocity dispersion in these more massive galaxies 
would make their spectroscopic analysis of weak absorption 
lines impossible. We have further trimmed their sample of 
137 galaxies by rejecting the ones which had the lowest S/N 
and the emission line galaxies, leaving us with a total of 101 
galaxies. 

Table |l| shows the distribution of Andreon et al. (1996 ) 
and Andreon et al. (1997) morphological types for the C93 
galaxies. Immediately obvious is the fact that 12 of the Cald- 
well early-type galax ies are in fact late-types according to 
Andreon et al. ( |1996; ). Although only one of these turns out 
to lie more than two standard deviations blueward of the 
U - V CMR ridge line (see Fig. [j] and §3.1). It must be 
noted, however, that up to date morphological information 
for many of the fainter galaxies in the sample is not avail- 
able, and that this is where most of the blue galaxies reside. 
Nevertheless, only three out of the seven bluest galaxies have 
no determined morphologies. 



t An ASCII table of the galaxy colours and line indices is avail- 
able on request via email from A. I. Terlevich. 
t We take the re cessional velocity of th e Coma cluster to be 
~ 6800-fCms - 1 (Colless h Dunn 1996), and assume -ffo = 
WQKms~ 1 Mpc~ L , which gives a distance modulus of (m — M) = 
34.2. 

§ C93 converted the GMP (b-r) colours to fES-V) colours using 
galaxies in common with the compilation of ( Burstein et al. 1987|) 
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Table 1. The morphological typ es for the C93 galaxies (Andreon 



Figure 1. The U—V colour-magnitude relation for the subsample 
of 101 Coma galaxies for which we have appropriate S/N spectra, 
both colours and magnitudes are taken using 8.8 arcsec diame- 
ter apertures. The standard deviation about the best fit (solid 
line) is 0.07 ± 0.01 mag. The symbols are coded by each galaxy's 
morphological type. Galaxies which lie between 1 and 2 standard 
deviations blueward of the CMR best fit have an additional circle 
drawn around the graph marker. Galaxies which lie more than 2 
standard deviations blueward of the CMR best fit have two ex- 
tra circles drawn about their graph markers. The same symbol 
scheme is used for the galaxies in all the figures throughout this 
paper. 
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Figure 2. The index— index relation between the two metal line 
indices used in this paper for the sample of Coma galaxies. A 
grid of single age stellar population models from Worthey (1994) 
is overplotted. Solid lines connect points of equal age; dashed lines 
connect points of equal metal abundance. The figure shows that 
for absorption strengths of C24668 > 6 A, the models do not 
match the parameter space occupied by the observations very 
well. This is most simply interpreted as an 'overabundance' of 
C24Q68 compared to Fe. The symbols, which are defined in Fig- 
ure pi are coded by galaxy morphology, and residual from the 
CMR. 



Andreon et al. 199' 



), split into the three coarse bins 



of Spiral/Late type, SO and Elliptical (see text). 



Morphological type 


Number of galaxies 


No morphology 


31 


Spiral/Late-type 


12 


SO 


47 


E 


11 


All types 


101 



2.2 Spectroscopy 

The C93 spectra were all taken on the KPNO 4m, using the 
HYDRA multi-fibre positioner. They have a spectral resolu- 
tion of 3.8 A (FWHM), and the fibre aperture was 2 arcsec 
diameter. In order to co mpare the spec t ra with Worthey 



stellar pop ulation models (Worthey 1994; Worthey & Otta- 



viani 1997), the spectra were re-sampled to the Lick/IDS 
resolution by convolution with a Gaussian of wavelength 
dependant width. Because the C93 data does not include 
any observations of stars in common with t he Lic k/IDS, the 
prescription used was that of Kuntschner (1998) who used 
it to correct his 4.1 A (FWHM) resolution Forna x data to 
the Lick/IDS system. Both the Kuntschner (1998) and the 
C93 spectra cover a similar wavelength range, and unlike 
the Lick/IDS spectra, whose spectral resolution degrades 
notably towards the blue (see Worthey & Ottaviani 1997, 
fig. 7) they have a constant spectral resolution over their 
entire wavelength range. The errors on the indices were esti- 
mated from the photon noise of the spectra, and do not take 
any systematic errors into account. The most likely source 
of systematic error is the re-sampling of the data onto the 
Lick/IDS system. Nevertheless, the possible systematic er- 
rors will have little effect on the relative comparisons made 
in this paper. We note that in addition, our Poisson error 
evaluation probably underestimates the true error because 
response variations from fibre to fibre will introduce extra 
scatter. 

The final step necessary to place the line indices into the 
Lick/IDS system is to correct the galaxy spectra for veloc- 
ity dispersion. Velocity dispersion correctio ns for individual 
line indices were taken from Kuntschner (1998), who cal- 
culated them by artificiall y broa dening stellar spectra (e.g. 
Davies, Sadler & Peletier (1993)). The C93 spectra are not 
of sufficient resolution and S /N to measure accurate central 
velocity dispersions for our galaxies, so instead we construct 
a 'Faber-Jackson' relation (Faber & Jackson 1976) between 
our V band 8.8 arcsec diameter aperture magnitudes my 



(see se ction 2.3), and the velocit y di spersions (a) of Lucey 
et al. (1997). A biweight fit (cf., §2.3) gives a relation of 



/0310(a) = (5.6 ± 0.6) - (0.22 ± 0.01) x m v 



(1) 



(with la bootstrap errors). We used this relation to find 
values of a for those galaxies in the C 93 sample which were 
not present in the Lucey et al. (1997) data. 



We use the H ^a and H8a Balmer line indices ( Worthey 
&i Ottaviani 1997 ) as the predominantly age sensitive spec- 
tral features. While Hf3 is more sensitive to age, it is also 
more affected by nebular emission, which can rapidly fill the 
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absorption feature (Gonzalez 1993), so is not used. Higher 



order Balmer lines are less sensitive to emission from ionised 



gas (Osterbrock 1989), making an accurate measurement of 
the true stellar absorption easier. In order to increase the 
signal to noise, we have used H8a + Hja in the final analy- 
sis. The C24668 feature is identified by Worthey (1994) as a 
particularly sensitive metallicity feature. However it is pos- 
sible that it suffers from over abu ndance problems in the 
larger, higher metallicity systems (Kuntschner 1998). Fig- 
ure [] shows the relation between C24668 and Fe4383 for the 
Coma galaxies, overplotted are the theoretical relation from 
Worthey (1994) models. Above a value of ~ 6 A, the mea- 
sured C24668 index values deviate from the region pre dicted 
by the models. A similar effect is seen by Kuntschner (1996 ) 
for early-type galaxies in the Fornax cluster. However in our 
data, the Fe4383 index has poorer signal to noise than the 
C24668 index, so we will use both, noting the discrepancy 
for the high metallicity galaxies. 



2.3 The photometry 

It is important to get as good a match as possible between 
the photometric and spectroscopic apertures, so that both 
are measuring the same part of the galaxy. In the case, for 
instance, of localised bursts of star formation occurring in 
the galaxy nucleus, the large photometric aperture and the 
small spectroscopic one, could well both be measuring dif- 
ferent mixes of stellar populations. 

In order to match the 2 arcsec diameter of the HYDRA 
fibres used for the spectra as best as possible, we use the 
8.8 arcsec diameter photometry from T99. The use of smaller 
metric apertures was not possible due to the increased ef- 
fects of variations in seeing conditions during the observa- 
tions. The T99 catalogue covers an area of 3360 arcmin 2 to 
a depth of V13 = 20 mag (ie., V band magnitude within 
a 13" diameter aperture) giving Johnson U- and V-band 
magnitudes for ~ 1400 extended objects. The RMS internal 
scatter in the photometry is 0.014 mag in V13, and 0.026 mag 
in U13 for V13 < 17 mag. 

Throughout this paper we use the biweight minimis- 
ing techni que, d escribed in T99 (see also Beers, Flynn & 
Gebhardt (1990)) to fit linear relations to the data. It was 
used for its resistance to outlying data points, and its robust- 
ness against non-Gaussian noise distributions. We derive the 
residuals of indices and colours from a mean index or colour 
- magnitude relation by the following equation: 



A(X) = X - (b x x m v + ex) 



(2) 



where X is a line index (C 2 4886, Fe4383 or HS A + Hf A ), or 
(U — V) colour, bx and cx are the slope and intercept of the 
best fit relation (see table ^|), and mv is the V magnitude. 



3 IMPLICATIONS FOR STAR FORMATION 
HISTORIES AND THE ORIGIN OF THE 
CMR 

3.1 The colour magnitude relation 

Given the selection criterion of C93, we should expect the 
CMR to be prominent in our (V, U — V) data-set. As can be 
seen in Figure [j], the bi-weight scatter is indeed small: a — 




C 2 466B (A) 

Figure 3. The distribution of galaxies in the (C24668, + 
Hta ) plane, wi t h a model grid of single stel lar population models 
( Worthey 1994 Worthey & Ottaviani 1997 ). Solid lines trace loci 
of constant age (1.5, 2, 3, 5, 8, 12 and 17 Gyrs from positive to 
negative H8 A + H^a values). Dashed lines trace loci of constant 
[Fe/H] ( -0.5, -0.25, 0.0, 0.25 and 0.5 from lowest to highest 
metal line strength). The different symbols correspond to each 
galaxy's morphology, and its offset from the colour magnitude 
relation as defined in figure [l]. 




Fe4383 (A) 

Figure 4. The distribution of galaxies in the (Fe4383, H<5a + 
H^va ) plane, wi t h a model grid of single stel lar population models 
( Worthey 1994 Worthey fe Ottaviani 1997 ). Solid lines trace loci 
of constant age, from 1.5 to 17 Gyrs. Dashed lines trace loci of 
constant [Fe/H] from —0.5 to 0.5. The different symbols corre- 
spond to each galaxy's morphology, and its offset from the colour 
magnitude relation as defined in figure pi. 



0.07±0.01 mag. Somewhat more surprising is the significant 
number of galaxies deviating blueward of the mean relation 
towards the faint end. The effect is not seen in the (B, B — 
V) relation of C93 (see figure 1 in C93), presumably due 
to the lower sensitivity of B — R colour to age, combined 
with their larger photometric errors. In order to trace the 
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position of these relatively faint and blue galaxies in the 
subsequent figures, we adopt a common labelling convention 
throughout the paper. The symb ols reflect the morph ologie s 
taken from Andreon et al. ( 1996 ) and Andreon et al. (1997), 
as defined by the key in figure [lj. The number or rings plotted 
around each symbol indicates its deviation from the best fit 
CMR. Points without rings lie either red-ward of the colour- 
magnitude relation, or are within la of it. Points with only 
one ring lie between la and 2a blueward of the relation, 
and points with two rings lie more than 2a blueward of the 
relation. Unless otherwise stated, when we refer to 'blue' 
galaxies we are referring to the objects plotted with at least 
one ring, and 'red' galaxies are the objects plotted with no 
rings. 



3.2 Line Index Diagrams 

We have used the Fe4383, C 2 4668 and H5 A + H-y A indices 
to place the galaxies on Worthey's age/metallicity diagnos- 
tic diagrams (Figures |^ and ^). In these figures, solid lines 
trace loci of constant age from 1 5 to 1 7 Hyrs (positive to 



negativ e H8a + H~/a values'). Dashed lines trace loci of con- 
stant [Fe/H] from —2.0 to +0.5 (lowest to highest metal 
line strength). Note that small systematic offsets (Ss 0.5 A) 
between the model and the observed index systems may well 
be present. We must stress that these models represent sin- 
gle stellar population (SSP) systems, whereas our galaxies 
contain a mixture of stellar populations of different ages and 
metallicities. The indices we measure are therefore a lumi- 
nosity weighted mean index of all of the different popula- 
tions present in each galaxy, and the higher luminosity of 
a young stellar population can mean that even the addi- 
tion of a small percentage (by mass) of young stars into an 
old galaxy, can severely affect it's overall 'age'. For exam- 
ple, if one adds a 10% (mass), lGyr, [Fe/H] = +0.5 pop- 
ulation, to an 18Gyr, [Fe/H] = population, the resultant 
stellar population would have the same li ne strengths as 
a 2Gyr, [Fe/H] = +0.25 SSP ( [Trager 199^) . In this paper, 
when we refer to a galaxy as 'young', we are referring only 
to its luminosity weighted mean age. 

As we would expect, the blue and red galaxies occupy 
different areas of the line index diagrams. The blue galax- 
ies which deviate by more than 2a from the CM relation 
(Fig. |l|) tend to populate the low age portion of the grid. 
The trend is common to the lower colour deviation galax- 
ies as well: only six out of the 20 blue galaxies lie on the 
same part of the Worthey grid populated by the red galax- 
ies. Both figures imply a young age for the blue galaxies 
(an effect driven mainly by the H S A + H"/ A index) ; however 
they seem to indicate different trends for the the red galax- 
ies which make up the ridge of the CMR. In Figure @ red 
galaxies span a range both in metallicity and in age, seeming 
to indicate a far from simple formation scenario, especially 
as the more massive galaxies appear to have the younger 
ages. We must remember however, that the model grids at 
these high metallicities are far from secure. In particular, we 
have already noted the overabundance problem in the mod- 
els for galaxies with C 2 4668 > 6 Ain §2.2. If we ignore these 
galaxies, the span in age and metallicity of the red galaxies 
is reduced. 

Fe4383 does not suffer from the same overabundance 
effects as C24668, however it also has more of an age de- 



pendence, producing slightly more degenerate model grids. 
In addition, the measurements have larger relative errors. 
In figure ^, we can see again that the blue galaxies have 
younger ages than do the red ones, however they no longer 
exclusively populate the low metallicity portion of the grid. 
Some of the galaxies are in parts of the plot completely 
outside the parameter space covered by the model grid, al- 
though these are the galaxies with the largest error bars in 
the Fe4383 data. The red galaxies have a lower spread in 
metallicity, and do not follow the same trend in age as in 
figure |^, however, due to the increased scatter in the Fe4383 
data, and the increased degeneracy in the model grid, it is 
still not possible to say that they follow a single age popula- 
tion model line. This may simply be due to the lower signal 
to noise of the Fe4383 index. 

It is worth noting that there is a morphological segre- 
gation in figures ^| and ||. In both figures all of the elliptical 
galaxies populate the old portion of the model grids, while 
the SOs and late type galaxies are spread over the grid. This 
seems to indicate there are both old a nd young SO and late 
type galaxies, but only old ellipticals ([Kuntschner fc Davie: ' 



1998; ; [Kuntschner 1998| ; |Mehlert et al. 1998j ). The overall ef- 



fect seems to be connected with galaxy luminosity, with the 
lower luminosity galaxies having more varied star formation 
histories. 



3.3 Colour residuals vs line strength residuals 



In section 3.2 



we showed that the blue galaxies which devi- 
ate from the mean CMR have younger luminosity weighted 
ages than do the red galaxies. In this section we shall further 
investigate the correlation between colour offsets from the 
CMR, and analogous offsets from spectral index vs. magni- 
tude relations (IMR) (Figure [|. 

We define the mean IMR in the same way as we defined 
the mean CMR, by using the biweight estimator (see sec- 
tion 2. J) to get a linear relation between the spectral line 
strength and luminosity of a galaxy. We list the results of 
the fits in Table ^, and they are shown as the solid lines in 
figure ^[ The residuals from the mean IMR are defined using 
equation [| 

Figure ^| and Table [| show that we can define IMRs 
for all of our spectral indices. The C2 4668 relation (Fig. |E|a) 
shows the best correlation, which fits in with the interpre- 
tation of the CMR as a metallicity sequence, since C24668 
i s the most met allicity sensitive of our three spectral indices 
( |Worthey 199^ ). The H8 A + H-y A relation (Fig. |c) also has 
an obvious IMR. Unlike the C2 4668 IMR, the red galaxies 
make up a relation of increasing H5 A + Hy A with decreasing 
luminosity, while most of the blue galaxies are deviant from 
the relation. 

The Fe4383 IMR (Fig |b), like the H8 A + H-y A IMR, 
shows increased scatter from the blue galaxies. To demon- 
strate the effectiveness of the biweight minimisation fit, we 
show some examples of different types of fit to the data. 
The solid line is the biweight scatter minimisation fit to all 
the data. The dashed line is the biweight fit to only the red 
CMR galaxies, the dot-dash line shows the biweight fit to 
galaxies with my < 16, and the dotted line shows the fit to 
all the data, using a method which bisects the ordinary least 
s quares fits made by minimising the X a nd the Y residuals 
(|lsobe et al. 1990|; |Feigelson fc Babu 1992j). The biweight fit, 
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Table 2. Results of regression analysis with la bootstrap errors for the line indices and colour vs V-band magnitude (see figures 
and 



Index 


Intercept 


Slope 


Scatter 


U-V 


3.5 ± 0.7 mag 


-0.14 ±0.01 


0.08 mag 


C 2 4668 


28.9 ±3.7 A 


-1.45 ±0.06 kmag- 1 


1.05 A 


Fe4383 


7.4 ±8.6 A 


-0.14 ±0.07 kmag' 1 


0.70 A 


H5 A + Hj A 


-18.4 ±9.3 A 


0.75 ± 0.12 kmag' 1 


1.44 A 



and the least squares fit to only the red galaxies are almost 
identical. As is the case with the HSa + Hja relation, most 
of the deviation comes from the blue galaxies, so we again 
take the biweight fit as the definition of the mean popula- 
tion. It should be noted that there is a much larger bootstrap 
uncertainty in the slope of the Fe4383 IMR than with any of 
the others (see Table Q), making the existence of a non-zero 
slope to the relation only a 2a result. 

Out of the three indices we investigate, C24668 is the 
most sensitive to metallicity, and indeed shows the greatest 
correlation with galaxy luminosity. This diagram agrees with 
the standard paradigm of a metal abundance driven CMR. 
The poorer correlation of Fe4383 is initially surprising, but 
seems to reflect the lower metal abundance sensitivity of 
the Fe4383 index relative to C24668, a difference that is en- 
hanced by the so-called 'over abundance' effect (cf., Fig. |^). 
This leads to the interesting possibility that the CMR is in 
fact driven by a correlation between luminosity and metal 
overabundance, however further investigation of this possi- 
bility would require improved spectra, and is beyond the 
scope of the present work. 

Finally, we note that models of the Fe4383 index have 
more of a dependence on age than C24668; thus an age gra- 
dient along the CMR (in the sense of fainter galaxies having 
younger ages) would tend to steepen the slope. This effect 
is seen only in the position of the blue galaxies in Figure [|b. 

We now investigate how deviations from the CMR cor- 
relate with deviations from the IMRs. We show the correla- 
tions in Figure ^| To aid the eye, we show the ±lcr deviations 
from the CMR and IMRs as dashed horizontal and vertical 
lines. The dotted line in panel (c) shows the ±3cr sca tter i n 
the HSa + Hja relatio n. We have used the Worthey ( 1994 ) , 
Worthey & Ottaviani (1997) models to show how changes 
in age (solid arrow) and metallicity (dotted arrow) would 
affect galaxies in the CMR/IMR plot. The age arrow was 
calculated by adding a 15% (by mass) population of lGyr 
old stars with solar metallicity, to an 8 Gyr old population 
with solar metallicity. The metallicity arrow was calculated 
by changing the metallicity of an 8 Gyr old population from 
[Fe/H] = to [Fe/H] = -0.5. These arrows can be thought 
of as vectors, in that increasing the size of the post-starburst 
population, or decreasing the metallicity of the population 
further, tends to increase the length of the arrows, without 
significantly changing their direction. Note however that if 
a galaxy were undergoing a starburst at the time of obser- 
vation, it would have Balmer lines (Fig^c) in emission, and 
would therefore have negative A(H8a + Hja), but it would 
still posses blue colours (negative A(U — V)). This emission 
quickly fades and turns into absorption, and after ~ lGyr, 
the galaxies indeed roughly follow the vector shown in Fig- 
ure pp. 



Figure ga shows deviations from the C24668 IMR. Only 
three of the A(U — V) < —2a (double-ringed) galaxies 
have AC24668 values which lie outside the ±la range indi- 
cated by the dashed lines, with the other four double-ringed 
galaxies lying within ±l<r of the AC24668 distribution. The 
age (solid) vector, indicates that we wouldn't expect much 
C 2 4668 deviation for A(U - V) = -0.3, and indeed the 
majority of the double-ringed blue galaxies lie close to the 
age vector. The three double-ringed galaxies most deviant 
in C24668, do lie closest to the metallicity vector, however 
given the large range in the AC24668 indices of the 'red' 
population, they are not inconsistent with the age vector. 
By contrast, the age and metallicity vectors for the AFe4383 
relation (Figure pb) are degenerate, and indeed all but a 
couple of the blue galaxies lie along these vectors. 

The correlation between A(U — V) and A(H5a + H^a) 
is the strongest of the three (Fig [i^) . Like Figure []a, there 
is a significant separation between the age and metallicity 
vectors, and all but one of the double-ringed blue galaxies 
lie closer to the age vector than the metallicity one. This 
includes the three blue galaxies which lie along the metal- 
licity vector in Figure 0a. In this case however, the range in 
the A(HSa + H^a) indices for the 'red' galaxies is less than 
the equivalent range of AC24668 in Figure |^a, and cannot be 
used to account for the distance of most of the 'blue' galaxies 
from the metallicity vector. Additionally, all of the double- 
ringed blue galaxies deviate from the mean HSa + H^a IMR 
by more than 3a (the dotted horizontal line). These effects 
can only be accounted for by using a younger mean stellar 
population. 



4 IMPLICATIONS 



We have shown using both the Fe4383 and the C24668 in- 
dices, that the 'red' galaxies which make up the colour- 
magnitude relation in the Coma cluster, span a range of 
metallicities. Whether this is the sole driving force behind 
the CMR, or whether there is also a correlation between age 
and luminosity is more debatable. There is no evidence for 
an age trend with the age-estimates based on the Fe4383 
index, however the C24668 models indicate that the more 
luminous galaxies have younger ages. It could be that the 
effect is not visible in the Fe4383 models due to the higher 
degree of degeneracy in the age and metallicity tracks for 
this index, but it should be noted that any trend in age 
comes from the C24668 index, at strengths where overabun- 
dance effects mean that it is no longer properly predic ted b y 
the Worthey models (see F ig [| and Vazdekis et al. (1996), 
(|L997| )). Kuntschner ( 1998| ) corrects the C 2 4668 index in his 
Fornax cluster galaxies for such overabundance problems in 
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Figure 5. The relation between the line indices and V magni- 
tude. The solid line in each panel shows the biweight regression 
fit used to calculate A(X) (see Table ||). The different symbols 
correspond to each galaxy's morphology, and its offset from the 
colour magnitude relation (see Figure hi). In panel b, the dashed 
line is the biweight scatter minimisation fit to only the red galax- 
ies, the dot-dash line shows the biweight scatter minimisation fit 
to galaxies with ray < 16, and the dotted line shows the fit to 
all the data, using an ordinary least squares bisector method (see 
main text). 
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Figure 6. The deviation of the galaxies from the mean relation 
in the three spectral indices (A(C24668) in panel a, A(Fe4383) in 
panel b, A(H8a + ^7a) m panel c), against their deviation from 
the mean relation in colour (A(C7 — V)). The fit to the mean rela- 
tions are calculated from Figures |lj and ^| and are summarised in 
table ^| The different symbols correspond to each galaxy's mor- 
phology^ and its offset from the colour-magnitude relation (see 
Figure hh . The dashed horizontal and vertical lines represent the 
±lo" dispersion in the colours and line strengths. The dotted hor- 
izontal line in panel c shows the 3<r dispersion in H8a + ^7a- 
The two arrows show how a change in galaxy age (solid arrow) 
and a change in galaxy metallicity (dotted arrow) would affect 
a 'normal' galaxy. The age arrow is calculated by adding a 15% 
population of 1 Gyr old stars to an 8 Gyr old galaxy of the same 
(solar) metallicity. The metallicity arrow shows how an 8 Gyr old 
solar metallicity galaxy would move if its metallicity were changed 
to [Fe/H] = -0.5. 
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the models, and shows that much of the age upturn for the 
bright ellipticals disappears. If the correlation between age 
and luminosity is real however, it would be in the direction 
predicted by semi-analytic hierarchical clustering models, 
which predict that more luminous cluster galaxies should 
have younger luminosity weighed ages, and higher metal- 



liciti es, than the less luminous ones (Kauffmann & Chariot 
^c^rnally, we note that the sense of this age gradient is 
such that the CMR cannot be explained by an age variation 
alone. In contrast, even greater metal abundance variations 
are required in order to restore the observed slope. 



We find that galaxies which deviate bluewards from 
the CMR, also deviate from the (HSa + #7a) — magnitude 
and the Fe4383-magnitude relations, but not the C24668- 
magnitude relation. This is consistent with the colour de- 
viant galaxies having undergone a secondary burst of star 
formation giving young luminosity weighted ages, or show- 
ing overall younger ages. The interpretation of the CMR as 
a mainly metallicity sequence, with deviations being caused 
by age variations seems correct. This allows the CMR to be 
used to place constraints on the star formation histories of 
t he cluster galaxies , an approached used by many authors 
(Bower et al. 1992| ; Stanford et al. 1998^ ~ 
Bower et al. 1998| ). 



Ellis et al. 1997 



Despite the attempts of C93 to reject late-type galaxies 
from their sample, a few still seem to have slipped in. We 
found that the SO and late-type galaxies span a range in 
ages (|Kuntschner fc Davies 1998|; |Mehlert et al. 199g|), but all 



elliptical galaxies are old. In addition, it is important to note 
that the younger galaxies are also lower luminosity galaxies, 
an effec t also seen in the Fornax sample of Kuntschner & 
Davies (|l99cj), and in low lumin osity (Mb > 17. 5) Coma SO 
galaxies (jCaldwell fc Rose 1998;). Worthey (|l997j) also found 



that lower luminosity galaxies have a larger spread of ages. 
This may indicate that star formation activity is being 'down 
sized' to lower luminosi ty objects as the universe becomes 
older (Cowie et al. 1996) or may simply be an artifact of the 
greater numbers of lower luminosity galaxies. Of the galaxies 
which lie on the CMR (the red galaxies) we find no obvious 
difference between the ages of the early and late types, all of 
them have old luminosity-weighted mean stellar ages. This 
includes 10 out of the 12 late-type galaxies in the sample, 
although this result should not be too surprising since their 
spectroscopic properties will almost certainly be dominated 
by their bulge. 

Finally, we note that the use of CMR-corrected (U — V) 
colours is an efficient way of searching for post-starburst 
galaxies. The (B, B-V) CMR constructed by C93 for these 
same objects (figure 8 in C93), shows that them to have a 
much smaller deviation from the mean population than our 
U—V colours do. In fact all our very blue objects, which de- 
viate blueward by more than two standard deviations from 
the CMR colour, have Balmer line strengths that deviate 
by more than three standard deviations from the red galax- 
ies. If we were to search for the PSB galaxies from scratch, 
and made a colour cut one standard deviation blueward of 
the CMR, we would find all of the 3<r deviant Balmer line 
galaxies, with only a 50% contamination from the 'normal' 
population. 
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